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ABSTRACT 



Aims. We want to understand the chemistry and physics of disks on the basis of a large unbiased and statistically relevant grid of disk 
models. One of the main goals is to explore the diagnostic power of various gas emission lines and line ratios for deriving main disk 
parameters such as the gas mass. 

Methods. We explore the results of the DENT grid (Disk Evolution with Neat Theory) that consists of 300 000 disk models with 1 1 
free parameters. Through a statistical analysis, we search for correlations and trends in an effort to find tools for disk diagnostic. 
Results. All calculated quantities like species masses, temperatures, continuum and line fluxes differ by several orders of magnitude 
across the entire parameter space. The broad distribution of these quantities as a function of input parameters shows the limitation 
of using a prototype T Tauri or Herbig Ae/Be disk model. The statistical analysis of the DENT grid shows that CO gas is rarely 
the dominant carbon reservoir in disks. Models with large inner radii (10 times the dust condensation radius) and/or shallow surface 
density gradients lack massive gas phase water reservoirs. Also, 60% of the disks have gas temperatures averaged over the oxygen 
mass in the range between 15 and 70 K; the average gas temperatures for CO and O differ by less than a factor two. Studying the 
observational diagnostics, the [C n] 158 j-im fine structure line flux is very sensitive to the stellar UV flux and presence of a UV excess 
and it traces the outer disk radius (i? ut). In the submm, the CO low J rotational lines also trace R oM . Low [O i] 63/145 line ratios (< 
a few) can be explained with cool atomic O gas in the uppermost surface layers leading to self-absorption in the 63 fim line; this 
occurs mostly for massive non-flaring, settled disk models without UV excess. A combination of the [O i] 63 line and low J CO lines 
correlates with several disk properties such as the average O i gas temperature in disks, the outer disk radius, and the UV excess. 
Conclusions. The [Oi]63/CO 2-1 line ratio is a powerful diagnostic to break disk modeling degeneracies. A combination of the 
[Oi] 63 yum flux and the [Oi] 63/CO 2-1 line ratio can be used for M gas < 10~ 3 M Q to obtain an order of magnitude estimate for the 
disk gas mass purely from gas observations. The previously used conversion of a CO submm line flux alone generally leads to larger 
uncertainties. 

Key words. Astrochemistry; circumstellar matter; stars: formation; Radiative transfer; Methods: numerical; line: formation 



1. Introduction 

In the recent decade, the modeling of the gas chem- 
istry and physics in protoplanetary disks has become a 
field of its own, following the lead of the dust SED 
modeling with multi-dimensional co n tinuum radiative trans- 
fer (e.g. IChiang & Goldreichl Il997b iD'Alessio etail Il998t 
iDullemond et alJ l2002h . The impact of the gas chemistry 
and energy balance is largest in the disk surface layers 
where most of the infrared and sub-mm line emission is 
arising (e.g. iKamp & Dullemondl 12004 IJonkheid et alJ I200H: 
iGorti & Hollenbachll2004 iNomura & Millar! 120051) . The com- 
plexity of the gas chemistry and individual heating/cooling 
processes as well as their mutual dependence makes the gas 
structure modeling more complicated, physically more demand- 
ing, and computationally more time consuming. Hence, most 
modeling has so far concentrated on either single selected ob- 
jects (e.g. lOi et al.ll2003t I Semenov et al 1 120051; lOi etail 120081; 
iGorti & Hollenbachl2008t Kenning et al.ll2010h prototypes (e.g. 



Aikawa et all 120021; iMeiierink etail 120081; iGoicoechea et alJ 
20091) or small model series studying the dependencies of 



specific parameters on t he gas emission from disks (e.g. 
l Aikawa & Nomura! 120061: IJonkheid et all 120071; INomura et all 
l2007HGlassgold et al.ll2009UKamp et alJl2010h . 

Current and upcoming observing facilities such as Herschel, 
ALMA, JWST/MIRI, SPICA/SAFARI are excellent tools to 
study the gas in large samples of protoplanetary disks. Hence, 
there is a clear need for a simultaneous and consistent modelling 
of dust and gas in protoplanetary disks, with the aim to iden- 
tify the diagnostic power of certain continuum and line observa- 
tions for the physical, chemical and temperature structure in the 
disks. These identifications must be based on a large variety of 
disk models covering a wide spread in stellar, disk, dust, and gas 
parameters. 

We jointly produced a grid of 300000 disk models ((Disc 
Evolution with Neat Theory) to study the predictive power 
of individual gas lines, line ratios and continuum tracers 
dWoitke et al.l l2010). The line emission studies concentrate on a 
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series of fine structure (O i and C n) and molecular lines (high 
J CO and H2O) that are observable with Herschel as well as 
a number of CO sub-mm lines accessible from the ground. 
A first analysis of the grid has shown that for example the 
[0 1] 63 /im emission strong ly depends on the a mount of UV ex- 
cess coming from the star (Woitke et ai1 l2010l) . Subsequent ap- 
plication to the fir st findings from th e Herschel Open Time Key 
Program GASPS dPinte et al .11201 Oh revealed that the measured 
[0 1] fluxes for a sample of 30 disks indicate that disks around 
T Tauri stars indeed require extra heating of the disk surface — 
either through a UV excess and/or X-ray irradiation — , while 
the emission from disks around Herbig Ae stars can generally be 
matched by photospheric heating alone. This paper now presents 
a more systematic and complete study of the gas properties and 
diagnostics from the DENT grid. The dust diagnostics will be 
discussed in a separate paper. 

In the following, we briefly su mmarize the DENT g rid ap- 
proach (Section [2]i as presented in I Woitke et al.l (l2010h before 
we discuss the statistics of gas chemistry and thermal structure of 
the 300 000 models in Sect.[3]and the line diagnostics in Sect.|4] 
We present our current best strategies for deriving disk gas pa- 
rameters in the final conclusions (Sect. |6j. 



2. The DENT grid 

The DENT grid h as been described for the first time in 
Woitke et al. (2010) and we summarize here only a few key as- 
pects. More details can be found in the appendix. 

The grid is computed us ing the 3D Mon t e Car lo radia- 
tive transfer code MCFOST (iPinte etalj [20061 l2009h and the 
gas chemistry and physic s modules of ProDiMo dWoifke et alJ 
l2009at iKamp et al]l2010h . The column densities, flaring, dust 
size distribution and settling are parametrized with power laws, 
assuming a prescribed global gas-to-dust ratio. Mie theory is 
used to compute the dust optical properties assuming an astro- 
nomical silicate composition. The central star is parametrized 
by its mass and age, as well as strength of the UV excess. Dust 
and gas temperatures are calculated separately from Monte Carlo 
continuum radiative transfer and a detailed heating/cooling net- 
work, respectively (see Sect. l2.1l and l2.2) . 

We have chosen 1 1 free model parameters for the grid, lead- 
ing to a total of 322 030 disk models. The range of parameter 
values explored in DENT is discussed further in the appendix. 
The total computing time of the grid is 200 000 CPU hours, an 
average of 40 CPU minutes per model. The total grid took 3 
weeks on 400 processors and generated 1.8 TB of data. 

The grid is very coarse in nature, e.g. stepping through gas 
mass in factors of ten. Hence, we cannot expect any method that 
we develop to measure disk properties such as the disk gas mass 
from gas line emission or line ratios to be more accurate than that 
step size, in this case an order of magnitude. The same holds for 
other quantities such as the outer disk radius, the flaring index, 
the dust-to-gas mass ratio etc.. 

For the following analysis, it is also important to keep in 
mind that the grid of disk models is not constructed to reflect 
the statistics with which disks occur in nature, e.g. the fact that 
some objects are rarer than others. There could be some disk 
types that we missed in our approach, or others that simply do 
not occur in nature. We take the latter into account by exclud- 
ing certain parameter combinations in the analysis of the grid. 
Models with a flaring index of /3 — 0.8 are for example more 
appropriate for late stages of disk evolution such as debris disks 
and so they are excluded from the study of line diagnostics in 



younger disks (noted in the figure captions whenever appropri- 
ate). A global disk analysis using the Toomre criterium shows 
that 94% of the disk models in the Dent grid are gravitationally 
stable (see Appendix lA.3t . 

We chose parametrized density structures to avoid im- 
posing the complex gas structure that is based on presump- 
tions mai nly derived from disk mo d eling approaches such as 
th ose of iGorti & Hollenbachl (120041) . iNomura & Millarl (120051) 
or IWoitke et alJ l2009a). While it is true that the more con- 
sistent vertical disk structure models need less assumptions 
than the simple parametrized models used in DENT, even 
a full ProDiMo in cluding the vertical structure calculation 
dWoifke et al.l l2009al) is not guaranteed to capture all relevant 
physics. Using parametrized disk structures enables us to ex- 
plore a large unbiased parameter space. From these simplified 
disk structures, we can still assess the relative influence of sev- 
eral of the key parameters on line and continuum fluxes. Even 
though the modeled fluxes are not absolute, we can still test fun- 
damental aspects of our disk understanding, such as uncertain- 
ties in various observational techniques to derive disk dust and 
gas masses. 



2.1. MCFOST 

MCFOST is 3D Monte Carlo continuum and line radiative trans- 
fer codaJHie detail s of the numerical sch emes are presented 
in iPinte et all d2006h and iPinte etafl d2009h . In short, the code 
stochastically propagates photon packets in 3D through the disk. 
The transport of packets is governed by successive scattering, 
absorption and re-emission events, which are determined by 
the local dust properties (opacity, albedo, scattering matrix) and 
temperature. MCFOST u ses the immediate re-emission concept 
(Biorkman & Wood 2001) with a c ontinuous d eposition of en- 
ergy to estimate the mean intensity (lLucvll99 9). In the optically 
thick regions of the disk, the code uses a diffusion approximation 
method to converge the temperature structure. Emerging SEDs 
are generated by a ray-tracing run from the temperature and ra- 
diation field estimated by the Monte Carlo run. 

For the DENT grid, 10 6 packets were used to compute the 
temperature structure. The radiation field was converged using 
10 4 packets per wavelength bin at wavelength larger than 0.5 /-im. 
As a reliable estimate of the UV field is critical for the chemistry 
module in ProDiMo, we used 10 5 packets per wavelength bin 
for wavelengths shorter than 0.5 jjm. This ensures that the UV 
radiation field is well converged, even in the deep regions of the 
disk. 

MCFOST also includes a NLTE line transfer module which 
is not fully used here, but whose ray-tracing module is used in- 
stead to compute the emerging spectral lines from the level pop- 
ulations estimated by ProDiMo. 



2.2. ProDiMo 

ProDiMo calculates radiation thermo-chemi cal models of proto- 
pla netary disks. It is de scribed in detail in Woitke et al. (2009a) 
and lKamp et alJ d2010l) and we summarize here the most salient 
features. Most noticable, the code is not used in its full capacity, 
because the self-consistent computation of the vertical hydro- 
static equilibrium structure is omitted. 

The code takes as input the dust density and tempera- 
ture, nd ust (r, z), Td ust (r,z) and dust properties as provided by 
MCFOST. In addition, the mean radiation field J v (r,z) from 
MCFOST is passed on to ProDiMo. There are two grid param- 
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eters that apply only to the gas part of the model: the gas mass 
M gas , and hence the dust-to-gas mass ratio 6, and the fractional 
UV luminosity defined as fyy = Luv /L* in the wavelength range 
between 91 and 250 nm. 

The chemical network consists of 9 elements, 71 species 
connected through 950 reactions (neutral-neutral, ion-molecule, 
photoreactions, cosmic ray reactions and adsorption & desorp- 
tion of CO, C0 2 , H 2 0, NH 3 , CH 4 ). The background radiation 
field J v enters the computation of the gas chemical and energy 
balance at several points, namely in the detailed frequency de- 
pendent integration of the UV ionization and dissociation cross 
sections, in the radiative pumping rates involved in the non-LTE 
modeling of atoms and molecules (O i, C i, C n, Mg n, Fe n, Si n, 
S ii, CO rotational & ro-vibrational, 0-H2 & P-H2 ro-vibrational, 
0-H2O & P-H2O rotational), and in the photoelectric heating 
rates. The atomic and molecular data used in the statistical 
equilibrium calculat ion is described in detail in Sect. 6.1.5 of 
IWoitke et al. (2009a|). 



3. Gas physics and chemistry 

As a first step in the statistical analysis of the DENT grid, we dis- 
cuss here how the masses and average temperatures of a number 
of key species (O, CO, C + and H2O) depend on the grid parame- 
ters (see Appendix lA.4t . With this method, we can obtain a first 
understanding of the disk chemistry (atomic versus molecular), 
the relevance of ices (how much CO or H2O is bound in ice 
form) and the thermal disk structure (gas dust de-coupling and 
cool surface layers). 



3.1. Oxygen and CO masses 

Oxygen is the third most abundant element in the interstellar 
medium. In dense molecular clouds or protoplanetary disks, oxy- 
gen can be either in the gas or ice phase (both are included in our 
chemical model); in the gas phase, it can be atomic (O) or molec- 
ular (CO, O2, H2O). By measuring all possible phases, we can 
constrain the total oxygen mass in a disk. 



Table 1. Maximum possible fraction of total disk mass for se- 
lected species X sp assuming all carbon or oxygen are locked into 
that species. Also given are the species mass (m sp ), abundance 
(e) and the expected values from the statistical distribution of 
gas mass fractions. 




-4.0 -3.5 -3.0 

og Mass(O) [Msun] / Mgos [Msun] 

Fig. 1. The atomic oxygen mass fraction X(0) = M(0)/M gas for 
all disk models with M gas > 10~ 7 M . Color coded is the total 
disk gas mass. Note that the y-axis is logarithmic. 
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Fig. 2. The CO gas mass fraction X{CO) = M(CO)/M gas for all 
disk models with M gas > 10~ 7 M . Color coded is the total disk 
gas mass. 
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Fig. Q] shows the statistical distribution of atomic oxygen 
masses. Low mass disk models have almost no CO, H2O or ices 
(X(CO) < 10 3 , X(U 2 0) < 10" 4 ) and hence their atomic oxy- 
gen mass fraction is close to maximum. Assuming all oxygen is 



atomic, we obtain a mass fraction of X(0) = 2.3 x 10 3 . Higher 
mass disk models show a wider distribution of atomic oxygen 
mass fractions in Fig.[T] This is due to a higher molecular frac- 
tion (dust shielding and gas self-shielding), e.g. CO, CH4, and 
H2O. Studying the CO mass fraction distribution (Fig. |2), the 
number of disks with masses larger than 10~ 5 M Q that have a 
particular CO mass fraction is rising all the way towards the 
maximum possible CO mass fraction. The lowest mass disks 
(10 -6 < M gas < 10~ 7 M ) show instead a broader and flatter 
distribution; ice formation and photodissociation for those mod- 
els are extremely model dependent. Also, none of those models 
ever shows all carbon being bound in CO. The dust-to-gas mass 
ratio lowers for example the peak of the CO mass fraction by a 
factor 5 going from 0.001 to 1 (Fig. [3). This is due to the average 
dust temperature being lower in the high dust-to-gas models and 
hence more CO freezing out. 
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Fig. 3. The CO gas mass fraction X(CO) = M(CO)/M gas for all 
disk models. Color coded is the dust-to-gas mass ratio. 



Fig. 5. Mean gas temperature of oxygen (see Eq.tfTJi) for all mod- 
els. Color coded is the stellar luminosity. 
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Fig. 4. The H 2 mass fraction X(H 2 0) = M(H 2 0)/M gas for disk 
models with an inner radius equal the sublimation radius. Color 
coded is the surface density power law index e. 



3.2. Water masses 

The mass fraction of gas phase water depends strongly on the 
surface density power law index, i.e. on the distribution of mass 
within the disk. Massive disk models that concentrate the mass 
in the inner dense disk (e = 1.5), convert a substantial fraction 
of their oxygen into water (Fig.|4]i. For shallower density gradi- 
ents, the outer disk contains more mass than the inner disk; here 
temperatures are low so that water is predominantly frozen on 
the grain surfaces. What remains is a low fraction of gas phase 
water above the ice layer. Models that have an inner gap, i.e. 
Ri n = 10 or 100 times the sublimation radius tend to have lower 
mass fractions of gas phase w ater. Those models l ack the big 
inner gas phase water reservoir dWoifke et al.ll2009bh . 



3.3. Oxygen gas temperatures 

We calculate average gas temperatures for various species from 
the models following 



,(r,z)TJr,z)dV 



(1) 



where n sp is the species density and m sp its mass. M sp is the total 
mass of a species in the disk model. The oxygen gas tempera- 
ture distribution peaks at 20 K and the median gas temperature 
is 40 K. The distribution shows that ~ 60% of the models stay 
between 15 and 70 K. The temperature distributions for specific 
subsamples of parameters, such as the the stellar luminosity and 
disk gas mass, are shown in Fig. [5]and[6] If we split the sample 
by luminosity, the higher luminosity models show on average 
higher gas temperatures. This can be directly attributed to the 
higher energy deposition rate which leads to an increased gas 
heating. Looking at low mass disk models (Fig. [6]), we note that 
their average gas temperature is almost a factor three higher than 
that of high mass disk models. This is explained by optical depth 
effects: in models with moderate optical depth, a larger fraction 
of the total mass can be heated by irradiation. Hence, models 
with lower disk mass tend to have on average higher gas temper- 
atures. These general trends are based on a huge number of disk 
models and the actual distribution of gas temperatures is fairly 
wide with a long shallow tail towards very high gas temperatures 
(Fig.©. 

The O i gas temperatu re statistics (and emerging line fluxes) 
have already been used in lPinte et al.l ( 2010l) to identify the need 
for an extra heating mechanism from the [Oi] line fluxes of 
T Tauri stars. In that particular case, the stellar luminosity of 
the low mass stars does not provide enough energy to explain 
the observed flux levels. Our models show that a UV excess 
(/uv = 0.1), which can originate either from a chromosphere 
or from ongoing accretion, can provide this extra energy and re- 
produce the observed fluxes. X-ray heating is an alternative pro- 
cess for T Tauri disks (|Nomura et al1l2007 ; iGorti & Hollen bach 
120081: lAresuet al.ll2011l) . 

Fig. [7] shows the impact of dust settling on the coupling of 
gas and dust temperature for the regions in which atomic oxygen 
dominates. Since atomic oxygen is widely distributed in the disk, 
we can take this to first order as an estimate of the amount of 
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Fig. 6. Mean gas temperature of oxygen (see Eq.([T)) for all mod- 
els. Color coded is the disk gas mass in three bins representing 
optically thin disk models, disk models in the transition from 
thin to thick, and optically thick disk models. 
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Fig. 7. Ratio of gas and dust temperatures of neutral oxygen as a 
function of the dust settling parameter (s = - without settling, 
s = 0.5 - with settling). This ratio is a measure of the temperature 
de-coupling in the disk surface. 



thermal gas-dust de-coupling. Well mixed models, 5 = 0, show 
a fairly narrow peak around T g = Tj and a tendency to have gas 
temperatures higher than dust temperatures (often in the form of 
a hot surface gas layer). The settled models, s = 0.5, show on av- 
erage gas temperatures that are lower than those of the dust. This 
is not due to the presence of cool surface layers in settled models 
(see Sect I3.6t ; the mass averaged Oi gas temperature is hardly 
affected and so the reason lies in the average dust temperature. 
In settled models, only the small grains stay behind in the disk 
surface, while the larger ones settle towards the midplane (see 
Table IA.ll . Hence, the effective grain size in the surface of set- 
tled models is smaller compared to non-settled models, thereby 
decreasing the dust emissivity at long wavelength. This leads to 
slightly warmer dust. 
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Fig. 8. Mean gas temperature of CO (see Eq.([T)) for all models. 
Color coded is the outer radius R oa . 
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Fig. 9. Mean gas temperature of CO (see Eq.([T)) for all models. 
Color coded is the surface density power law exponent e. 



3.4. CO gas temperatures 

Fig.[8]and|9]show the distribution of CO mass averaged gas tem- 
peratures in the disk models as a function of the outer radius and 
the surface density power law index. Given our choice of surface 
density power laws, it is evident that the CO mass is dominated 
by the material close to the outer radius. Since the gas tempera- 
ture is also a clear function of distance, models with larger outer 
radius show on average lower CO gas temperatures. Fig. |8]also 
illustrates the limitations of the concept of a single canonical 
CO temperature Tco to convert submm line fluxes into disk gas 
masses. We will get back to this point when looking for possibil- 
ities to estimate disk gas masses from line fluxes in Sect. l5.2l and 
I5.3I The surface density power law exponent also impacts the 
distribution of mass in the disk. Hence, models with a low e put 
more mass in the outer disk compared to the inner disk; hence, 
they show on average lower CO gas temperatures. In addition, 
the trends with luminosity and disk mass found for oxygen, hold 
also for CO. 



6 



I, Kamp et al.: Continuum and line modeling of disks around young stars 




log T 9 (0!)/T g (CO) 



Fig. 10. Difference between average gas temperatures of atomic 
oxygen and CO for all models. Color coded is the dust-to-gas 
mass ratio. 



For a large fraction of models, the average gas temperatures 
for atomic oxygen and CO are very similar. Except for the lowest 
mass disks, both temperatures are within a factor two of each 
other. This also holds across all gas-to-dust mass ratios as can be 
seen from Fig. [10] 



3.5. Ionized carbon in disks 

Ionized carbon forms a thin skin on the disk surface with a typi- 
cal column density of 10 I7 ~ 18 cm~ 2 . Very low mass disk models, 
M gas < 1CT 5 M G , are so optically thin that carbon is fully ionized 
throughout the disk resulting in a constant M(C + )/M gas fraction 
of ~ 1.2 x 10~ 3 . This is only weakly dependent on detailed disk 
parameters. 

As soon as the ionizing UV radiation is blocked by suffi- 
cient dust opacity, carbon turns atomic and subsequently molec- 
ular. The higher mass models show thus a decreasing fractional 
C + mass, with the total M(C + ) leveling off at ~ 1(T 8 M Q . The 
variation between models can be an order of magnitude mostly 
depending on the choice of disk extension and stellar luminos- 
ity (Fig. ITTb : the C + mass is higher for larger outer disk radii 
and higher luminosity stars. The gas residing in the outer disk is 
generally more optically thin — also to interstellar UV radiation 
— than the gas in the inner disk regions. 

The dependency of the C + gas temperatures on stellar lu- 
minosity is less pronounced than for atomic oxygen. The distri- 
bution has a strong peak at low temperatures around 20 K ex- 
cept for luminosities L„ > 3.6 L Q (Fig. \l7\i . There is a clear 
dependency on the outer radius similar to CO, with large models 
(^out = 300, 500 AU) having an average C + gas temperature of 
<7^ + > ~ 25 K. The small models (R ont = 100 AU) show a much 
broader flat distribution between 25 and 100 K. The reason is 
that the C + mass averaged temperature is dominated by the out- 
ermost regions — just as for CO — since it is evenly distributed 
over the entire disk surface. Hence, the inner radius plays no sig- 
nificant role here. 
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Fig. 11. Mass of ionized carbon as a function of total disk gas 
mass. Color coded is the outer disk radius. 




Tg_CI 

Fig. 12. Mean gas temperature of ionized carbon (see Eq.(Q3) for 
all models. Color coded is the stellar luminosity. 

3.6. Cool surface layers 

A subset of DENT grid models shows very cool surface layers 
on top of warmer midplanes, i.e. models with a positive gas tem- 
perature gradient (as seen from the surface). These are very flat, 
massive disks in which the grains have started to settle. Typical 
models in this category show M gas > 0.01 M , no UV excess, 
settled dust (s = 0.5) and are non-flaring (J3 - 0.8, 1.0). These 
models appear here for the first time because of the wide unbi- 
ased parameter study. 

A related effe c t has been seen earlier in models by 
lAikawa & Nomural (120061) who studied grain growth in disk 
models that do solve iteratively for the vertical hydrostatic equi- 
librium. For well mixed models with a large maximum grain size 
flmax, those authors find a similar effect. Photoelectric heating is 
diminished in those models resulting in lower gas temperatures 
in the disk atmosphere. On the other hand, optical depth is re- 
duced as well allowing stellar UV photons to penetrate deeper 
and heat the gas and dust at intermediate heights. The authors 
also find that disks with larger a max flare less than those with 
small molecular cloud like grain size distributions. This is an 
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Fig. 13. Statistical distribution of [Cn] 158 fim emission as a 
function of disk gas mass. Color codes is the logarithmic disk 
gas mass. 



Fig. 14. [Cn] 158 fjm fluxes as a function of stellar luminosity 
for disk models with M gas > 10~ 5 M . Color coded is the UV 
excess /uv- 



immediate result of the lower gas temperatures in the surface 
layers. For their largest maximum grain size a max = 10 cm, 
they find that the gas in the disk atm osphere becomes c ooler 
than the dust. In a subsequent paper iNomura et al.l (120071) in- 
cluded the effect of grain settling and X-ray heating and show 
that — compared to the standard well-mixed model with molec- 
ular cloud like grains — this indeed lowers the gas temperature 
in the surface and increases the gas and dust temperatures close 
to the midplane. The gas temperature gradient becomes much 
flatter, but there is no pronounced positive temperature gradient. 
Unfortunately we cannot disentangle from their published model 
the effect of additional gas heating by X-rays and decrease in gas 
ph otoelectric heating due to dust growth and settling. 

iJonkheid et al.l y007) have studied a series of four Herbig 
disks with an increasing gas-to-dust mass ratio due to grain 
growth and settling. The basi c dust disk models are taken from 
iDullemond & Dom inik ( 2004), who use t he dust midplane tem- 
perature to set the scale height of the disk. lJonkheid et all (12007) 
do not iterate the vertical structure based on the gas tempera- 
ture and so their approach is somewhat similar to our choice of 
parametrized disk models. They find that a large fraction of the 
disk gas mass is indeed cooler than the dust and that the highest 
gas-to-dust mass ratios lead indeed to a positive gradient in the 
vertical gas temperature. 

In nature, we do not expect such flat, massive disks, with 
strong settling to be very common. However, the presence of 
such cool surface layers would have implications for the disk 
surface chemical composition and line emission. At cold tem- 
peratures, neutral-neutral chemistry becomes inefficient and the 
formation of warm water in the disk surface could be affected. 
Self-absorption in the dominant cooling line [O i] 63 /mi (see 
Sect. 14.3b and temperature gradients as measured from various 
12 CO and I3 CO lines could reveal such cool surfaces. 

4. Line diagnostics 

In the following sections, we explore the diagnostic power of 
individual emission lines from the DENT grid. We check espe- 
cially for clean dependencies on disk input parameters such as 
outer radius, disk mass, dust-to-gas mass ratio, surface density 
power law exponent etc. 



4. 1. The C + fine structure line at 158 fim 

Since ionized carbon originates from a thin skin around the 
disks, its mass (column densities) tends to be constant until the 
disk gas mass itself becomes so low that the disk is transparent 
to UV photons. The total [Cn] 158 pm flux correlates well with 
the mass in C + indicating that the line is mostly optically thin. 

Fig.Q~3]illustrates this effect nicely through the statistical dis- 
tribution of [C n] 158 yum fluxes as a function of disk mass. Down 
to disk masses of 10~ 4 M , the disks stay optically thick and the 
C + mass fraction and hence its emission stays constant around 
~ 3 10~ 18 W/m 2 . The large width of the flux distribution for these 
models reflects the spread in stellar luminosity, UV excess, outer 
disk radius, flaring angle; the higher fluxes results from models 
with high luminosities (stellar and/or UV — see Fig. [T4l . large 
outer radii, and flaring geometry. The [C n] line flux depends 
only weakly on the disk surface density gradient, the gas-to-dust 
mass ratio, settling, the minimum grain radius a m j n and the inner 
disk radius. For lower gas masses, the peak of the [C n] distri- 
bution shifts towards lower fluxes and the distribution becomes 
much narrower. 

The DENT grid was generated using only two values of UV 
excess, 10% and 0. 1 % of the total stellar luminosity. While these 
values bracket the UV excess observed for young T Tauri stars, it 
can lead to an overestimate of the UV radiation for Herbig Ae/Be 
stars that do not possess an excess in the first place. Fig. [14] 
shows that the [C n] emission increases by an order of magnitude 
for the high UV models. The Herbig Ae/Be models with pure 
photospheric UV fluxes would fall even below the red points in 
Fig.[H 

In practice, the [Cn] 158 /urn line is easily contaminated 
by the low density surrounding gas, even with Herschel/PACS. 
Although the line is potentially a good tracer as shown above, it 
requires a very careful observational procedure and data reduc- 
tion to ensure that the flux measured is indeed purely from the 
disk. 

4.2. Oxygen fine structure lines at 63 and 145 fim 

In the DENT grid, 80% of the disk models show a 145 /im line 
that is a factor 10 - 100 weaker than the 63 //m line. Given that 
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the Herschel sensitivity limits for both lines are very similar, this 
cle arly favors the dete ction of the 63 /mi line. 

Woitk e et alJ {2010) have shown that the 63 //m fine struc- 
ture line could be promising as an order of magnitude disk gas 
mass estimator. Their Fig. 4 illustrates that if we can estimate the 
average excitation temperature of neutral oxygen, the measure- 
ment of the 63 /mi line can be used to derive gas masses in the 
context of a statistical analysis of a large sample of disks such as 
the GASPS sample. 

4.3. The [Oi] 63/145 line ratio 

In PDRs, low [Oi] line ratios (< 10) are expect ed if both lines 
are optically thick (iTielens & Hollenbachl [l985) and the gas is 
cooler than ~ 200 K. iLiseau et alJ d2006l) argue that confusion 
with cool envelope gas might explain rati os of a few in obser- 
vations of YSO's. lLorenzetti et alJ d2002l) suggested that high 
optical depths in the [Oi]63 and 145 /mi lines can explain the 
[Oi] 63/145 ratios of 2 - 10 observed by ISO/LWS in massive 
disks around early-type HerbigAeBe stars. 

The grid results presented here confirm that the [O i] 63/145 
ratios for massive (>10 2 M ) disks is of the order of 1 - 10. 
In our models, we observe that [Oi] 63/145 < 10 can be pro- 
duced by the disks themselves, when [O i] 63 /mi is about to go 
into absorption, whereas [Oi] 145 /mi is still in emission. As a 
fundamental line, [O i] 63 /mi can go into central absorption if 
the emitted light has to pass through another cold neutral oxy- 
gen layer above, whereas the [Oi] 145 /im line practically never 
goes into absorption. These effects happen for massive disks 
only, in particular for settled dust distributions, and are strongly 
dependent on inclination, see Fig. Q3] and also Sect. 13.61 The 
hypothesis that such models exist is a unique testcase for the 
SOFIA/GREAT instrument that can observe the [O i] 63 //m line 
with high spectral resolution and thus check for self-absorption 
within the line profile. 

The median [Oi] line ratio for the entire DENT grid is 25. 
Contrary to PDRs, the [Oi] 63/145 line ratio of disk models is 
not sensitive to the average oxygen gas temperature over a wide 
range of temperatures between 50 and 500 K. The line ratio cor- 
relates instead with e.g. the dust-to-gas mass ratio (Fig. |T6}. 
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Fig. 16. Statistical distribution of [Oi] 63/145 line ratios as a 
function of the dust-to-gas mass ratio. 



Fig. 17. [Oi]63 /mi versus 12 CO J=2-l line emission. Color 
coded is the outer disk radius with values of 100, 300 and 
500 AU. 



4.4. The CO low J rotational lines 

The low J 12 CO rotational lines are optically thick and arise 
mainly from the outer disk surfaces. Hence they form an ex- 
cellent tracer of the outer disk radius. The larger the disk, the 
larger the emitting surface. The oxygen fine structure line at 
63 /um is not affect ed by the disk outer radius for R out > 100 AU 
dKamp et alJl2010h . Fig. ITTl shows the mean [Oi] line flux plot- 
ted against the CO 2-1 flux with the outer disk radius coded in 
color. In this diagram, the disk models clearly separate into small 
ones with R out = 100 AU and ones that are larger than 300 AU. 
This indicates that a combination of the oxygen and CO line can 
help to remove the outer disk radius degeneracy in measuring 
disk masses. 

Detection limits with current radio telescopes such as APEX 
or SMA are a few times 10~ 20 W/m 2 . With ALMA coming 
along soon, we will be able to detect much fainter disks down 
to ~ 10~ 22 W/m 2 . Fig [18] shows the CO 2-1 line flux versus the 
dust continuum at 1 .2 mm with the respective ALMA band 6 de- 
tection limits. These grid results suggest that ALMA is sensitive 
enough to detect disk masses of the order of 10~ 6 M . 

4.5. The[Oi]63/C0 2-1 line ratio 

Fig. [T9l illustrates that the [Oi]63/CO 2-1 line ratio correlates 
with the average oxygen gas temperature in the range 15-70 K. If 
the [O i] 63 /mi line is optically thin, it depe nds on disk gas mass 
and average oxygen gas temperature (Woi tke et al.l 12010). The 
CO 2-1 line is always optically thick (down to M gas ~ 10~ 4 M ) 
and is therefore a function of the CO gas temperature and disk 
surface area, i.e. ^ ou t- Using the two lines [O i] 63 /mi and CO 2- 
1, we also get a handle of the amount of extra UV radiation irra- 
diating the disk (see Fig. [20b. 

From Fig. Q3 [19] and [20] it becomes clear that this [Oi]/CO 
line ratio captures in a complex way a number of disk properties 
such as the average gas temperature, the disk outer radius, the 
flaring index and the amount of extra UV irradiation (/uv)- Thus 
this line ratio is a powerful means to break the modeling degen- 
eracy and find a suitable method to estimate gas masses even in 
the absence of a large number of observing constraints. This is 
especially relevant for large surveys such as GASPS (Dent et al. 
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Fig. 15. Unusual behaviour of the [O i] lines in some disk models. For massive disks with settled dust distribution, the fundamental 
[Oi] 63 fim line may go into central absorption (l.h.s.) whereas the [Oi] 145 fim line stays in emission (r.h.s.), which complicates 
the data interpretation. Model parameters M* = 2M , age = 3Myr, (r eff = 4960 K, L* = 3.66 L ), /uv = 0.001, M gas = 0.1M G , 
dust/gas = 0.01 , R m = 17.8 AU, R out = 500 AU, e = 0.5, /3 = 1 .0, s = 0.5, a min = 1 fim. 
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Fig. 18. 12 CO J=2-l line emission versus continuum flux at 
1 .2 mm. For discrete intervals in CO line flux and 1 .2 mm flux, 
mean disk masses are calculated from the ensemble of grid mod- 
els within that bin. The colored contours show the distribution of 
these mean disk gas masses. Overplotted are the line and contin- 
uum sensitivity limits of ALMA band 6. 



2011 in preparation) that contains many objects where only the 
[O i] 63 fim fine structure line is detected and ancillary data is 
mostly restricted to SED's and CO submm lines. 



4.6. Water lines 

Due to the molecular structure, water emission lines are not 
only difficult to calculate (non-LTE, population inversion, opti- 
cal depth effects), but also difficult to interpret in terms of global 
disk parameters such as gas mass or surface density profiles 
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Fig. 19. Histogram of the [Oi]63/CO 2-1 line ratio in mod- 
els with 15 < <r g (0)> < 70 K (60% of all models, see 
Sect. [33}. Color coded is the mass average oxygen gas temper- 
ature <r g (0)>. 



(ICernicharo et al.ll2009l) . IWoitke et all (l2009bl) traced the origin 
of water lines in the disks around Herbig Ae stars back to three 
main water reservoirs: (1) a main water reservoir in the central 
midplane out to ~ 10 AU (containing most of the water mass), 
(2) a cold water belt at r w 10 - 100 AU and z/r ~ 0.05 (origi- 
nating from desorption of water ice mantles) and (3) a hot water 
surface layer at distances of 2 - 25 AU (originating from an ac- 
tive neutral-neutral chemistry at T > 200 K). Water lines that 
originate deeper in the disk suffer from the problem that gas and 
dust temperature couple at some fiducial depth, making it poten- 
tially impossib le to see the line emission above the stron g dust 
continuum dWoitke et alJl2009bt lPontoppida n^raIll20irj|) . 
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Fig. 20. [Oi]63 /mi versus 12 CO J=2-l line emission. Color 
coded is the amount of UV excess, /uv- 



Fig. |2TI shows that for massive disks (M gas > 0.01 M ), the 
water lines are generally stronger in the settled models. The 
effect of dust settling on the water chemistry and line emis- 
sion is very complex as it affects the available dust surface 
area, the temperature gradients and the continuum optical depth. 
Disentangling those effects simply from the grid observables is 
impossible. Hence, we discuss here two exemplary models to 
illustrate the general trend displayed in Fig. |2T] Picking a set- 
tled and non-settled flaring disk model with a m ; n = 1.0 /mi 
(M, = 2 M , T eS = 4960 K, 3.7 L , fl in , out = 17,500 AU, 
M g . ds = 0.1 Me, 8 = 0.01, e = 0.5,0 = 1.0), we find that the dust 
in the surface layers of the settled model is warmer on average, 
while the gas temperature barely changes (using EqJT]for oxy- 
gen, it changes from 13 to 16 K). This was already discussed in 
Sect EH 

Even though the temperatures are low, there is very little 
grain surface area available for e.g. H2 or ice formation. The 
H2 formation on grain surfaces is thus inefficient and atomic hy- 
drogen abundances are high in these settled models. As a conse- 
quence of the lack of grain surfaces in the settled model, molec- 
ular hydrogen forms via the H route and the atomic hydrogen 
abundance stays a factor 100 higher than in the non-settled mod- 
els. This leads to an efficient cold water formation route driven 
by radiative association 



H + O 
H + OH 



OH + hv 
H 2 Q + hv 



(2) 
(3) 



The relevance of radiative association for cool atom ic H regions 
was already pointed out bv lGlassgold et al] (2009), but found to 
be irrelevant for the inner disk structures that they studied. These 
reactions have low rate coefficient^], but with atomic hydrogen 
being a factor 100 more abundant in settled models, they domi- 
nate the formation route of water. The surface layers (z/r > 0.2) 
of the settled model show water abundances that are several or- 
ders of magnitude higher than those of the unsettled model. 



1 The UMIST rates are k(R + O) = 9.90 10^' 9 (r/300 K)" 
and £(H + OH) = 5.26 10- 18 (77300 K)~ 522 exp(90/T) 
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Fig. 21. Ortho-H20 179 /mi line flux versus continuum 160 /mi 
flux for massive disk models with M gas > 0.01 M . Color coded 
is the settling parameter: 5 = denotes homogeneously mixed 
models, s = 0.5 denotes models in which the dust grains have 
settled to a smaller scale height H(r, a) oc H(r) ( 005 m )~° 25. 



5. Diagnostic tools 

In the following, we discuss a few methods to measure the disk 
gas masses from gas line observations. It is evident by the way 
how we set up the parameter space that the DENT grid does not 
reflect a proper distribution of disk properties, e.g. number of T 
Tauri versus Herbig disks, fractional distribution of disk masses 
reflecting disk lifetimes etc.. Hence, our aim is not to decide 
which method is better, but to put the various methods into the 
context of disk modeling and discuss some of the limitations 
and uncertainties that become evident if one confronts the simple 
methods with a wide variety of disk appearances (sizes, flaring 
indices, gas-to-dust mass ratios, etc.). 

The most common technique to estimate disk gas masses has 
been the use of the mm dust continuum flux. At those wave- 
lengths, the disk is assumed to be optically thin. Hence, the 
mm-flux should reflect the entire dust mass contained in large 
grains. The main uncertainties in converting the mm-flux into 
a total dust and subsequently gas mass are the dust opacities 
k v ~ ko(v/vo)~P and the dust-to-gas mass ratio 5. An average 
opacity index /? ~ 0.46 has been rece ntly found for the disks 
in the p Ophiuchus star forming region riRicci et al.ll2010l) : these 
low values of (3 (compared to 1 .7 for the dust in the interstellar 
medi um) are gener ally interpreted in terms of grain growth, see 
also iDrainel (|2006). However, recent interferometric data sug- 
gests that p varies with location in the disk (e.g. Ilsella et all 
2010). The dust-to-gas mass ratio is generally assumed to be the 
canonical value of 6 — 0.01 found for the interstellar medium 
neglecting potential differences in dust and gas evolutionary 
timescales. Even if the dust mass could be estimated very pre- 
cisely, this method would still leave us will large uncertainties on 
M gas because of the unknown conversion factor. Hence, the in- 
trinsic uncertainties of this method are much larger than an order 
of magnitude. 

Very few direct estimates of the ga s-to-dust mass ratio exist 
for disks. An example is the study by dGlauser et al.ll2008t) that 
estimates dust and gas mass on the same disk line-of-sight in a 
disk. The derived ratio of 220^ [12 is compatible with the canon- 
ical value of 100. 
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Fig. 22. Two colour line flux diagram, [Oi] 63/145 versus 
[Oi]63/[Cii] 158, of a sub-selection of 10740 DENT models 
with pd/p g = 0.01, e = 0.5, s = 0.5, inclination angle 41°. Color 
coded is the disk gas mass. 



The following three sections discuss potential methods only 
based on gas emission line fluxes, their limitations and uncer- 
tainties. 

5.1. [Oi] and [Cn] fine structure line ratios 

Figure [22] shows a two-colour line flux diagram suggesting a 
method how to determine the disk gas ma ss purely from thre e 
line flux observations. As was argued by iKamp et"ail (|2010|) . 
the [Oi] 63/145 line ratio shows a clear correlation with gas 
mass. This plot includes the results for all types of central stars, 
low and high excess UV, flared and non-flared disks, continuous 
disks as well as disks with large inner holes, and varying dust pa- 
rameters and outer disk radii. Despite this variety of disk types, 
all models with a certain disk gas mass fall into certain regions 
in this two-colour diagram. Apparently, the main dependencies 
on /uv and j3 tend to cancel out here, because these parameters 
shift all line fluxes up and down equally. 

However, the diagram looks different for other sub-selections 
of models as the one stated in the caption of Fig. [22] A steeper 
surface density power law (higher e) weakens in particular the 
[Cn] 158 pm flux and hence shifts the points rightward in the 
figure. 

Selecting only settled models has shifted the location of mas- 
sive models (M gas > 0.01 M ) downward in Fig. 1221 compared 
to a similar plot for unsettled models. The reason for this is 
mainly [O i] 63 pm self-absorption in the cool surface layers as 
discussed in Sect. 13.61 

In addition, there are other parameters such as the outer disk 
radius R out and the dust size parameters, which introduce ad- 
ditional non-trivial dependencies and also a large scatter in the 
various mass bins. Thus a mass determination exclusively based 
on these three line ratios would implicate large errors especially 
for the higher mass disk models. 

5.2. Low J CO lines 

For many objects, the [Oi] 145 and [Cn] 158 pm lines will not 
be available as they are often at least a factor 10 fainter than the 



[O i] 63 pm line. Therefore, we study in the following also the 
CO low J lines that are accessible from the ground. Especially 
with ALMA coming up, these lines will be detectable for a large 
number of disks (see Sect l4.4l i. 

The low J 12 CO lines are often optically thick an d can hence 
provid e at most a lower limit to the total gas mass. iDent et all 
(2005) present a formula to convert the CO 3-2 line flux 7co3-2 
into a gas mass assuming LTE and a constant 12 CO/H2 conver- 
sion factor of 5 x 10~ 5 . We adopt here their approach and derive 
a relation between gas mass and intensity 

where 7co3-2 is the integrated CO 3-2 line intensity in K km/s. 
Since most of the CO mass will reside in the outer disk, the line 
flux will be dominated by cold gas. Hence, an average gas ex- 
citation temperature of r ex = 50 K is assumed in the following 
estimates. To apply Eq.© to our grid data, we need to convert 
the modeled integrated line fl uxes in W/m 2 to K km/s (formula 
adapted from lThi et al]|2004l) 

i3 

/co3-2 = 0.01 X Fcoa-a^fO- 1 . (5) 
2k 

Here, A = 8.6696 10~ 2 cm and k is the Boltzman constant in 
erg/K. The solid angle is assumed to be Q = 7r(HPBW/2) 2 , 
with HPBW = 13.7" the half power beam width of the JCMT at 
345 GHz. The formulas for the CO 2-1 line look similar; how- 
ever, the line is generally a factor 2.6 weaker than the CO 3-2 
line reflecting the difference in the Planck function (v3_2/v2-i) 3 . 

We insert our modeled fluxes first into Eq.(|5) and then the 
converted intensities into Eq.© to derive a gas mass estimate. 
These derived values underestimate the actual gas mass, espe- 
cially at the high mass end, where the CO 3-2 line is optically 
thick. Fig. [23] shows that the behavior is much more complex 
and that there is a regime, where the simple formula even over- 
estimates the disk gas mass, most likely due to the simplified 
assumption of a uniform excitation temperature of 50 K for all 
models. Most notably, the standard deviation across all disk pa- 
rameters for this method to measure disk gas mass is large, gen- 
erally more than two orders of magnitude. 

5.3. The [Oi] 63/CO 2-1 line ratio 

We proceed thus to suggest here an alternative approach using 
the insight on the diagnostic power of the [Oi] 63/CO 2-1 line 
ratio from the previous sections. Fig. [24] shows for each bin of 
[Oi] 63/CO 2-1 line ratios the resulting fit using a 3rd order poly- 
nomial. The coefficients for these polynomial fits are presented 
in Table [2] In the following, we check how accurate the disk 
masses derived from a combination of these two lines are. 

For every model, we calculate the [Oi] 63/CO 2-1 line ratio 
and then select the corresponding curve from Fig. |24] This curve 
in combination with the calculated [0 1] 63 pm flux results in a 
'measured' disk gas mass. This value is then compared to the 
actual gas mass in the models (Fig. [25j. The standard deviation 
indicates that this method can provide an order of magnitude 
estimate of the disk gas mass. The method does frequently over- 
estimate the disk gas mass by a factor two. Fig. [25] also shows 
that the relation levels off at disk masses beyond 10~ 3 M where 
the [Oi] 63 pm line becomes optically thick. 

For a large statistical sample in which individual disk pa- 
rameters are less well known, this method can indeed help to 
assess gas evolution in disks and to measure the gas dispersal 
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Fig. 23. Correlation between the gas mass estimated from the 
modeled CO 3-2 line fluxes (Eq.© and (0) and the original in- 
put gas mass of the models. The blue dots are the mean estimated 
gas masses and the error bar denotes the standard deviation. The 
red solid line is the expected correlation together with the dashed 
lines that denote an uncertainty of one order of magnitude. 
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Fig. 24. Correlation between the gas mass in the disk and the 
[Oi]63 /mi line. The second relevant parameter, the average 
gas temperature, is color coded through the [Oi]63/CO 2-1 
line ratio. A selection is made for CO 2-1 fluxes larger than 
10~ 23 W m -2 . The colored lines are polynomial fits of 3rd or- 
der to the data. 



timescale. The errors of this method are much smaller than the 
ones from the CO 3-2 line alone. Of course, more observational 
constraints on any object should naturally lead to more accurate 
measurements. 



6. Conclusions 

We have calculated a grid of 300 000 disk models using the dust 
Monte Carlo radiative transfer code MCFOST and the gas chem- 
istry and heating/cooling balance code ProDiMo. These models 
span an 1 1 -dimensional parameter space and enable us to study 
disk properties such as temperatures, species masses, gas line 



Table 2. Coefficients for the polynomial fit of M gas as a func- 
tion of [Oi] 63 /mi flux f^oi] and the [Oi] 63/CO 2-1 line ratio: 
A^gas = flo + ai-F[Oi] + a 2^ 0I| + fl 3^j 5 oi]- The fits are valid in the 
range -21 < F [0 i] < -14 W/m 2 . 
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Fig. 25. Correlation between the gas mass estimated (or 'mea- 
sured') from a combination of the [O i] 63 //m line flux and the 
[Oi] 63/CO 2-1 line ratio and the input gas mass of the disk. The 
blue dots are the mean 'measured' gas masses and the error bar 
denotes the standard deviation. The red solid line is the expected 
correlation together with the dashed lines that denote an uncer- 
tainty of one order of magnitude. 



fluxes and dust continuum fluxes in a statistical manner. The goal 
is to find suitable gas diagnostics that allow us to characterize the 
main disk parameters such as radial extent, gas mass and flaring 
from observational quantities. 

The main results found from an extensive statistical analysis 
of the entire DENT model grid are: 

- The [O i] 63 /mi traces disk mass but for massive non-flaring, 
settled disk models without UV excess, the line is optically 
thick and often suffers from self-absorption in upper cool 
disk layers. 

- The [Cn] 158 /mi fine structure line flux is very sensitive to 
the stellar UV flux and presence of a UV excess. 

- The CO low rotational lines trace the outer disk radius. 

- Dust settling enhances the water abundance in the surface 
layers through an efficient radiative association channel. The 
initial step of H2 formation proceeds in the gas phase via H . 

- Low [Oi] 63/145 line ratios (< a few) can be explained with 
cool gas in the uppermost surface layers instead of 63 /mi 
absorption by cool envelope remnant gas as suggested by 
iLiseau et al.1 (l2006h . 
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Table A.l. Parameters for the model grid and values assumed. 



stellar parameter 





stellar mass [M e ] 


0.5, 1.0, 1.5, 2.0, 2.5 


age 


age [Myr] 


1, 3, 10, 100 


L 


dependent parameter 


see Table [Ol 




dependent parameter 


see Table IA721 


/uv 


excess UV f V v = Luv/L* 


0.001, 0.1 


disk parameter 


M d 


disk dust mass [M Q ] 


10"', 10"°, 10"\ 10" 4 , lO" 3 




dust/gas mass ratio S 


0.001, 0.01, 0.1, 1, 10 


Rin 


inner disk radius [i? su bii] 


1, 10, 100 


Rout 


outer disk radius [AU] 


100, 300, 500 


e 


column density 






Nn(r) oc r~ e 


0.5, 1.0, 1.5 




flaring H(r) = H (±f 


0.8, 1.0, 1.2 


r 


reference radius [AU] 


100 


Ho 


scale height at ro [AU] 


10 


XlSM 


strength of incident ISM UV 


1 


<TcR 


cosmic ray H2 ionization rate 






[s- 1 ] 


5 x 10~ 17 


/pah 


abundance of PAHs 






relative to ISM 





a 


viscosity parameter 





dust parameter 


s 


settling 








0, 0.5 


"Win 


minimum grain size \jim\ 


0.05, 1 


^max 


maximum grain size \jim\ 


1000 


Pgr 


grain material density [g/cm 3 ] 


3.5 


radiative transfer parameter 


i 


inclination 


0°, 41.41°, 60°, 75.52°, 






90° (edge-on) 


v turb 


turbulent line width [km/s] 


0.15 



- The [Oi]63/CO 2-1 line ratio correlates with several disk 
properties such as the average 1 gas temperature in disks, 
the outer disk radius, and the UV excess. As such it is a pow- 
erful diagnostic to break disk modeling degeneracies. 

- A combination of the [0 1] 63 //m flux and the [0 1] 63/CO 2- 
1 line ratio can be used for M gas < 10~ 3 M to obtain an or- 
der of magnitude estimate for the disk gas mass purely from 
gas observations. The previously used conversion of a CO 
submm line flux alone generally leads to larger uncertainties. 
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Appendix A: The DENT grid 

The general outline of our numerical pipeline is described in 
IWoifke et all d2010h (see their Fig.l). In the first step, MCFOST 
is used to compute the dust temperature T<i ast (r, z) and frequency 
dependent radiation field J v (r, z) from radiative transfer. This in- 
formation is passed on in the second step to ProDiMo to compute 
the detailed gas chemistry n mo \{r,z) and temperature r gas (r, z). 
The resulting gas structure together with the non-LTE level pop- 
ulation numbers of the relevant species is inserted back into 
MCFOST to do the line radiative transfer and calculate the line 
profiles and integrated fluxes (step 3). In a second radiative trans- 
fer step, dust spectral energy distributions (SED) are calculated 
for all models (step 4). 
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Tabl e A.2. Stel l ar par ameters of the model grid interpolated 
from Siess et al for solar metallictities. 



M. [Mo] 


Age [Gyr] 


T eS [K] 


R, [Ro] 


L, [L ] 


l°gg 


0.50 


0.001 


3771 


2.13 


0.93 


3.5 


0.50 


0.003 


3758 


1.31 


0.36 


3.9 


0.50 


0.010 


3764 


0.85 


0.15 


4.3 


0.50 


0.100 


3828 


0.43 


0.05 


4.9 


1.00 


0.001 


4278 


2.62 


2.34 


3.6 


1.00 


0.003 


4262 


1.72 


1.00 


4.0 


1.00 


0.010 


4276 


1.18 


0.48 


4.3 


1.00 


0.100 


5357 


0.90 


0.66 


4.5 


1.50 


0.001 


4536 


3.11 


4.17 


3.6 


1.50 


0.003 


4600 


2.10 


2.00 


4.0 


1.50 


0.010 


5403 


1.84 


2.89 


4.1 


1.50 


0.100 


6872 


1.44 


4.63 


4.3 


2.00 


0.001 


4710 


3.60 


6.47 


3.6 


2.00 


0.003 


4961 


2.59 


4.12 


3.9 


2.00 


0.010 


8903 


1.69 


17.13 


4.3 


2.00 


0.100 


8829 


1.65 


15.94 


4.3 


2.50 


0.001 


4853 


4.13 


9.59 


3.6 


2.50 


0.003 


5697 


4.19 


18.61 


3.6 


2.50 


0.010 


10507 


1.81 


38.12 


4.3 


2.50 


0.100 


10217 


1.95 


39.72 


4.3 



A. 1 . Parametrized disk models 

For the shape and mass distribtion of the gas in the disk, we use 
the following parametric description 

between an inner and outer disk radius, R m and R out , respec- 
tively, with sharp edges. p(r,z) is the local gas mass density. 
The constant po is adjusted such that the integrated disk mass 
2n ff p(x, z) dz r dr equals Mdi s k- H(r) is the vertical scale height 
of the disk, assuming to vary with radius as 

H{r) = H [-\ . (A.2) 
\r / 

Ho is the reference scale height at reference radius ro. e is the 
column density powerlaw index and /3 the flaring power. 

The dust grains are assumed have a unique powerlaw size 
distribution 

dn(d) oc a ~ p da (A.3) 

between minimum grain radius a m ; n and maximum grain radius 
flmax- The free constant in Eq. (IA.3t is adjusted to result in the 
specifi ed dust/gas mass ratio pd/p (see Sect. 4.6 in Woitke et al. 
2009a), which is assumed to be constant throughout the model 
volume (e. g. no dust settling). The dust absorption and scatter- 
ing opacities are calculated by applying Mi e theory with optica l 
constants for astronomical silicate from iDraine & Leel dl984). 
The DENT grid comprises also models with dust settling, al- 
though that particular parameter is only sampled by the two ex- 
tremes of non-settled models and models in which grains larger 
than 0.05 /mi have settled to a smaller scale height than the gas, 
H{r, a) * H(r) (a/0.05 pm)"° 25 (for s = 0.5). 

Table lA. ll summarizes the free parameters in the DENT grid. 
For completeness, it also contains the fixed parameters ro, Hq, 

XlSM, %CR, /PAH, a, flmax, Pgr and Vturb- 

Effective temperature and luminosity are related parameters 
and are consistently c hosen from the ev olutionary tracks of pre- 
main sequence stars (ISiess et al . 2000) for solar metallicities. 
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Table lA.2l lists the parameter combinations chosen for the repre- 
sentative five stellar masses 0.5, 1.0, 1.5, 2.0, and 2.5 M Q . Disk 
models are calculated at four stellar evolutionary times, 1, 3, 
10 and 100 Myr. The stellar input sp ectra are interpolated from 
Phoenix stellar atmosphere models dBrott & Hauschildtl [2005 ) 
with solar metallicities. 

At UV wavelengths, however, where stellar activity and mass 
accretion create an UV excess with respect to the classical stellar 
atmosphere models, we switch to a powerlaw UV input spectrum 
with spectral intensity [erg/cm 2 /s/A/sr] 



Hoc A 



Puv 



(A.4) 



The flux is scaled to yield the prescribed fractional UV luminos- 
ity of / uv = Luv/i* = 0.1(high UV) or 0.001 (low UV) with 
the UV luminos ity Lyy being integ rated from 912 A to 2500 A 
as introduced by Woitke et al. (2010). Very few stars have a well 
known /uv (given the wavelength range defined above and the 
gaps in observational data of e.g. STIS and FUSE). One of the 
few stars is ET Cha, which has a value of 0.02 (Woitke et al. 
2011, submitted). 

A.2. The observables 

The quantities that can be compared to observations such as the 
gas line emission and dust SED are computed in step 3 and 4 
of the grid modeling procedure (see introduction to Sect.|2). All 
quantities are computed for 5 different inclinations that are uni- 
formly distributed, 0°, 41.41°, 60°, 75.52° and 90° (edge-on). 

The dust SED is computed at 57 wavelength points be- 
tween 0.1 and 3500 jjm. The line profiles and integrated fluxes 
of four species are computed: Cn, Oi, CO and H2O (ortho 
and para). Table IA.3I provides an overview of the respective 
atomic/molecular line data. 

A.3. Disk stability 

We choose here the Toomre criterium to check our disk models 
against gravitational instabilities. The Toomre Q parameter is 
defined as 

K(r)c s {r) 



Qkf) 



nGL(r) 



(A.5) 



where c s is the sound speed, "L the surface density, G the gravi- 
tational constant and the epicyclic frequency k is defined as 



K(r) 1 = 



1 d(r 4 Q(r) 2 ) 
r 3 dr 



(A.6) 



with Q. the orbital frequency. The disk is unstable at a certain 
distance r from the star if Q(r) < 1 . In the following, we derive 
an expression that provides a simple analytical way of checking 
'global' disk instability. For that, we make a few simplifying 
assumptions 

1 . The disks rotate with keplerian frequency fikep = -\/GM*/r 3 , 
hence /c(r) = £2 kep (r) 

2. The midplane gas temperature can be approximated by a ra- 
dial power law T = T^r/R^y 05 . 

We will discuss the impact of these simplifications further below. 

Inserting the power laws for temperature and surface density 
into Eq.( IA.5l ). we obtain 



Q(r) 



kT M t 



e-7/4 



Table A.3. Atomic and molecular data for lines calculated in 
the DENT grid. All data is taken from the LAMBDA database 
dSchoier et alj 120051) with references to the original data indi- 
cated in the last column. 



species 



Cn 



Oi 



o-IFO 



p-H 2 



A 



ident. 



reference 



157.74 



i pT i p 



[s" 1 ] 



2.300(-6) F77, L77, W02 



63.18 
145.53 



3 Pn- 3 P, 



8.865(-5) 
1.772(-5) 



B98, C80, J92, L77 



2600.76 

1300.40 

866.96 

650.25 

520.23 

433.56 

371.65 

325.23 

289.12 

260.24 

144.78 

90.16 

79.36 

72.84 



J=l-0 

J=2-l 

J=3-2 

J=4-3 

J=5-4 

J=6-5 

J=7-6 

J=8-7 

J=9-8 

J=10-9 

J=18-17 

J=29-28 

J=33-32 

J=36-35 



7.203(-8) 

6.91(-7) 

2.497(-6) 

6.12(-6) 

1.221(-5) 

2.137(-5) 

3.422(-5) 

5.134(-5) 

7.33(-5) 

1.006(-4) 

5.695(-4) 

2.126(-3) 

2.952(-3) 

3.638(-3) 



F01, J05, W06, Y10 



538.29 
180.49 
179.53 
174.63 
108.07 
78.74 



1 10 — > loi 

423 - > 3 12 
2l2 - * loi 
3o3 — * 2 12 
2l2 — » 1 10 
4 23 -> 3 12 



0.0035 

0.0306 

0.0559 

0.0505 

0.256 

0.484 



D02, F07, G93 



303.46 
269.27 
144.52 
138.53 
100.98 
89.99 



2o2 - > Fl 
In - > Ooo 
4i3 -> 3 22 

3 13 - > 2f)2 
220 — > Fl 
322 — > 2n 



0.0058 

0.0184 

0.0332 

0.125 

0.260 

0.352 



D02, F07, G93 



Refere nces: B98 dBell et al.1 [19981) . C80 dChambaud et alj [l980h . 
D02 dDubernet & GrosieaXl2002f). F07 dFaure et al.l 120071). F77 
dFlower&Launavl ll977Tr^F0 1 dFlowed 120011). G93 dGreen et al.l 
119931) , J92 dJaquet et alj 119921) . J05 dJankowski & Szalewiczl l 2005h . 
L77 dLaunav & Roueffl Tl977h. W06 dWernli et al.1 120061) . W02 
(Wilson & Bell 2002), Y10 dYang et alj2010h 



Here Tq and Xo denote the gas temperature and surface density 
at the inner radius of the disk Ri n . Since for our values of e (0.5, 
1 .0, 1 .5) Q(r) decreases with radius, disk models that obey Q > 1 
at the outer radius are also stable at smaller radii. Hence, in the 
following, we derive an expression for Q at r — R out , which can 
be used to determine the global stability of a disk model. 

The surface density and temperature at the inner radius can 
be derived from the total disk gas mass Mdi s k and t he luminosity 
and effective temperature of th e central star L„, T e ff (Tuthil let al.l 
l200lUvan der Plas et al.ll20TTh 



M dlsk (2 - e) 
AttR 2 



To 




Rout 
Rin 



(2-e) 



(A.8) 

(A.9) 



Using these, we can rewrite Eq. (IA.7t entirely in terms of the 
input parameters of the Dent grid 



(A.7) Q = 



8& 



e) 2 



R 



1/4 



(A. 10) 



I. Kamp et al.: Continuum and line modeling of disks around young stars 



15 



^out 



(2-e) 



^out 



(e-7/4) 



We apply this criterium to all 322 030 disk models and find 
that only 16422 of them are unstable according to the above 
derived criterium. Hence, 94% of all models are not affected 
by gravitational instabilities. This percentage changes by less 
than 1% if we (a) choose a more realistic temperature profile 
for optically thick disks T ~ (r / R m )~~°' 25 or (b) change To by 
a factor 2. The fact that Q(r) decreases with r still holds for 
T ~ (r/Ri n y°- 25 . Overall, we find that models with high gas 
disk masses — as expected — are more subject to instabilities. 
However, even in our highest mass bin the fraction of unstable 
models is only ~ 30%. The plots and conclusions of the anal- 
ysis of the DENT grid are not affected by the small number of 
models that could become unstable. 



A.4. Methodology 

The large number of disk models makes it impossible to carry 
out an individual study of their chemical and thermal structure. 
Hence in the following, we study the statistics of species gas 
temperatures and masses and how they depend on certain grid 
parameters such as the surface density distribution and settling. 

For the statistical analysis, we define a scalar quantity X, of 
any model i. It can be a parameter, or can be a measured value 
like the total mass of atomic oxygen M(O) of model ;. Next, 
we select an interval on which we perform the statistics [X^ - 
dx/2, X k + dx/2], where X k are sampling points for X. The total 
number of models (selected ;) in that selected range are N se i. The 
expected value and standard deviation of the distribution are then 
defined as 



(X) = J] Xi/N sel 



selected i 



cr 



V (Xi - (X)) 2 /N sel 

I selected i 



(A. 11) 



(A. 12) 
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